Once more about astrophysical S factor for the a + d — > 6 Li + 7 reaction. 
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Recently to study the radiative capture a + d — >• 6 Li + 7 process a new measurement of the 
6 Li(A 150 MeV) dissociation in the field of 208 Pb has been reported in [F. Hammache el al. Phys. 
Rev C82, 065803 (2010)]. However, the dominance of the nuclear breakup over the Coulomb one 
prevented from obtaining the information about the a + d — s> 6 Li + 7 process from the breakup 
data. The astrophysical 824(E) factor has been calculated within the a — d two-body potential 
model with potentials determined from the fits to the a — d elastic scattering phase shifts. However, 
the scattering phase shift itself doesn't provide a unique a — d bound state potential, which is the 
most crucial input when calculating the 5*24 (E) astrophysical factor at astrophysical energies. In 
this work we emphasize an important role of the asymptotic normalization coefficient (ANC) for 
6 Li — > a + d, which controls the overall normalization of the peripheral a + d —¥ 6 Li + 7 process and 
is determined by the adopted a — d bound state potential. We demonstrate that the ANC previously 
determined from the a — d elastic scattering s-wave phase shift in [Blokhintsev et. al Phys. Rev. C 
48, 2390 (1993)] gives 824(E), which is at low energies about 38% lower than the one reported in 
[F. Hammache et al. Phys. Rev C82, 065803 (2010)]. We recalculate also the reaction rates, which 
are also lower than those obtained in [F. Hammache et al. Phys. Rev C82, 065803 (2010)]. 
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PACS numbers: 26.35.+C, 25.45.-z, 25.40.Lw, 21.10.Jx 



The radiative capture a + d — > 6 Li + 7 is the only process that produces 6 Li in the big bang model. A special interest 
_L| I to this reaction has been trigerred by almost three order disagreemnet between the observational ratio 6 Li/ 7 Li [l[ and 
the calculated one [l[. Unfortunately direct measurements of the a + d — > 6 Li + 7 radiative capture are practically 
impossible at astrophysically relevant a — d relative kinetic energies E < 300 keV due to extremely low cross section. 
Hence, only indirect approach could be feasible to get information about 6 Li formation. The first indirect information 
about the astrophysical factor 824(E) for the a + d —> 6 Li + 7 process has been obtained in Q from the Coulomb 
breakup process 6 Pb( 6 Li (26A MeV), ad) 208 Pb. However, the energy behavior of the extracted astrophysical factor at 
low energies turned out to be constant what contradicted to all the calculations showing significant drop [3|. Recently 
in H a new attempt has been done to get the astrophysical factor 824(E) at astrophysically relevant energies from 
6 Pb( 6 Li (150AMeV), ad) 208 Pb. However, analysis has shown significant dominance of the nuclear breakup over the 
Coulomb one making impossible to determine the needed information about 824(E). Nevertheless, in Q the 824(E) 
has been calculated using the a — d two-body potential model. The potentials which are required to make such 
calculations were obtained from fitting the a — d elastic scattering phase shift for the s, p and d waves. The approach 
used to calculate the astrophysical factor is not related with the studied 6 Li breakup process. The only common 
information in the analysis of the breakup data and calculation of the astrophysical factor were the same bound state 
and scattering a — d potentials used to generate the corresponding bound state and scattering wave functions. In 
the potential approach used in Q the bound state potential, as we will discuss below, is the most crucial part of the 
input, which affects the overall normalization of the astrophysical factor. Unfortunately, the dominance of the nuclear 
breakup and dependence of the breakup data analysis on the optical potentials doesn't allow one to make a test of 
the quality of the adopted a — d potential. The approach applied in 0] is just repetition of the procedure used in 

Here we would like to discuss how reliable are the astrophysical factor 824(E) and the reaction rates for the 
a + d — > 6 Li + 7 radiative capture reported in Q and what should be done to improve our knowledge about them. 
It has been long ago recognized that the a + d — > 6 Li + 7 process at astrophysical energies is entirely peripheral 
reaction in the two-body potential model [3]. Evidently, the potential model itself is a limitation and it would be 
nice to check peripherality of this reaction within a many-body ab initio approach similat to what has been done 
recently for the 3 He + 4 He — > 7 Be + 7 in [6]. However, an important issue in a many-body approach remains to 
be solved is reproduction of the experimental binding energy, because the calculated astrophysical factor is sensitive 
to its value. Since such ab initio many-body calculations are not yet available we have to live with a more simple 
two-body potential model. The matrix element for the a + d — > 6 Li + 7 direct radiative capture in the long-wave 
approximation is given by 0] 

a _ a. I ^,r B Li 



M i f U, = ^ drl a % m) ^ , h (r). (1) 

All the notations are given in Q. The cut-off radius R c is introduced to reflect a peripheral character of the process. 
The multipolarity of the transition is A = 1, 2. It has been shown in [3J that the matrix element shows a very small 



sensitivity for R c < 4.5 fm, i.e. until distances which exceeds the 6 Li radius and makes possible to approximate the 
radial overlap function by the C a£ ;(oi) W— ij, 1/2(2 K a d,r)/r, where C^oi) is the asymptotic normalization coefficient 
(ANC) for the virtual decay 6 Li — > a + d in the state with relative orbital angular momentum I = and the 
total angular momentum J = 1, 1/2(2 n a( ir) is the Whittaker function determining the radial shape of the 

overlap function beyond of the a — d nuclear interaction region, r\ is the Coulomb a — d bound state parameter, 
K a d — a/2 Had^ad is the a — d bound state wave number, \i a d and e a d are the reduced mass and binding energy 
of a and d, correspondingly. As we can see, there are two inputs needed to calculate the matrix element. One is the 
a — d potential describing the continuum in the partial waves I = 1,2. This potential has been found in cite Q by 
fitting the elastic scattering phase shifts in these partial waves and it is a legitimate procedure. Note that at very low 
energies, say around the most effective energy of 70 keV, one can use a pure Coulomb scattering wave function in the 
initial state of the reaction. However, the role of nuclear interaction becomes very important with energy increase 
and it is responsible for reproduction of the resonance at Er — 0.7117 MeV. In the potential approach used in 0, H| 
the overlap function is replaced by the 

7 « d' (01) = S t!d (01) f<* d(ioi) (r) , (2) 

Here S a d (01) is the spectroscopic factor for the configuration (ac?)oi in 6 Li, <£ a d{wi){ r ) is the- radial wave function of 
the a — d relative motion in the field generated by the Woods-Saxon potential. Also n = 2 is the principal quantum 
number, i.e. the number of the nodes for r > is iV = 1, 2 = and J = 1 . The spectroscopic factor reflects the 
fact that the overlap function is not an eigenfunction of any Hamiltonian and, hence, is not normalized to unity in 
contrast to the bound state wave function. The potential, which is used to calculate the bound state wave function, 
has been determined in |4J from the fitting to the a — d elastic scattering phase shift in the channel I = 0, J = 1. Since 
the experimental elastic scattering phase shift includes the many-body effects of the scattered nuclei, the same is true 
for the two-body potential which fits the elastic scattering data. Hence, the spectroscopic factor in Eq. ([2]) should be 
set to S a d(oi) = 1- It is exactly what has been done in [4|. Since the reaction under consideration is peripheral at 
astrophysical energies E < 300 keV, functions in Eq. ^ can be replaced by their tales, i.e. 

r 6 Li r>R c r , W_ Vtl/2 (2K ad r) _ W_ ri , 1/2 {2n ad r) 

1 ad(01) ~ °Qd(01) — <W(ioi) , W 

i.e. in the two-body potential model in the asymptotic region C a d(oi) = b ad (wi), where 6 a d(ioi) is the single-particle 
ANC for the a — d bound state with the number of nodes for r > N — 1. The value of the single-particle ANC 
depends on the adopted bound-state potential. From the parameters of the bound state potential given in Q we 
find that b = 2.7 fin" 1 / 2 . Thus from the Woods-S axon potential given in Q we get the ANC, which is about 17% 
larger then the ANC used in Q. Since the astrophysical factor is proportional to the square of the ANC the usage 
of C 12 ^^ = 7.29 fm -1 rather than C^ d ( 01 ) = 5.29 fm -1 leads to the increase of the astrophysical factor compared to 
the one obtained in Q by almost 38%. It can be seen from Fig. 9 3, where both astrophysical factors are shown 
in the logarithmic scale. We came to the main point of this paper, the ANC, which is the most crucial input in the 
calculation of the 824(E). In [4| , as in [5[ , the s-wave scattering phase shift has been used to determine the bound-state 
potential, which generates the bound state wave function and, correspondingly, the ANC. However, it is well known 
from the inverse scattering problem that there is infinite number of the phase-equivalent potentials and to single out 
a unique potential one has to add two parameters (if only one bound state is present in the given partial wave): the 
binding energy and the ANC. In Q the adopted potential reproduces experimental a — d binding energy but still, 
one parameter, the ANC is still missing. Hence, the potential found in [4( is one out of an infinite set of the phase 
equivalent potentials and the ANC, which it generates, is not necessarily a correct one. However, there is another 
way of using the elastic scattering data which has been realized in Q ■ This approach is based on the analiticity of 
the elastic scattering S-matrix element what allows one to extrapolation the experimental scattering amplitude to the 
bound state pole in the momentum plane to get the residue, which is expressed in terms of the ANC. Once the ANC 
has been determined a unique potential can be found 0, Q ■ This potential has to satisfy the condition 

lim V N (r)e 2Kr -> 0, (4) 

r— *-oo 

where k is the bound state wave number. These procedures had been realized in 0] for the a — d s-wave elastic 
scattering. First, the Pade approximation was used to interpolate the elastic scattering amplitude in the physical 
region, which was extrapolated to the bound state pole to get the ANC. The obtained ANC C Q( i(oi) — 2.3±0.12 fm -1 / 2 
is lower than C^q!) = 2.7 fm -1 / 2 following from the potential used in After that the two-body potential has been 
found, which reproduces the s-wave a — d phase shift and provides correct binding energy and the same ANC. The 
found potential is more complicated than a Woods-Saxon one. It is a complex, energy independent potential, which 
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FIG. 1: (Color online) The red solid line is a — d bound state nuclear potential V (?) [4| and the blue dashed line is Vi^r). 
The explanation see in the text. 



can be expanded in terms of the harmonic oscillator basis. Evidently, that such a potential satisfies condition (0}. 
Thus the ambiguity in the two-body potential calls for more thorough selection of the potential because eventually 
the adopted potential for the bound state will determine the overall normalization of the astrophysical factor for 
peripheral direct radiative capture processes. 

If we would confine ourselves to astrophysical energies then it would be enough to replace in the matrix element 
(fTJ) the overlap function by its asymptotic term <j3j> and, evidently, the astrophysical factor will be proportional to 
^ad(oi)' However, if we want to extend our calculations to higher energies including the resonance region and above 
we need to carry more accurate calculations. To compare our results with the 824(E) factor presented in Q we adopt 
the same a — d scattering potential as in [ij . For the a — d bound state potential it would be logically reasonable to 
use the bound state potential from 0, H|. However, since it has quite complicated form, we simply modify the bound 
state potential used in Q using the theorem of the inverse problem in scatteringtheory [loj . This theorem allows one 
to recover a phase equivalent potential to the Woods-Saxon potential used in [4[ with arbitrary ANC. Assume that 
we find a nuclear potential V N (r), which together with the added Coulomb potential V c (r) fits the elastic scattering 
phase shift in the I = partial wave and the bound state wave function calculated with this potential has the ANC 
C . Then the phase-equivalent potential is given by 

V»{r) = V N (r) + **fL , (5) 



K(r) = Log[l + (r - 1)(1 - / dxx 2 ^p 2 (x))]. 

Jo 



(6) 



The bound state wave function in the potential V\(r) — Vj (r) + V c (r) can be expressed in terms of the bound state 
wave function (p(r) in the potential V(r) = V (r) + V c (r): 

^"" V M t^(,) ' (7 » 

From Eq. ([7]) taking the limit r->oowe obtain that C% = r _1//2 C. Let V(r) — V N (r) + Vc(?) be the nuclear Woods- 
Saxon plus Coulomb bound-state a — d potential adopted in Q , which generates the bound-state wave function with 
C a d(oi) = 2.7 fm -1 / 2 . Then for r = 1.378 we obtain the wave function in the potential V\(r), which has the same 
ANC as obtained in Q and used in In Fig. [T]both nuclear potentials are shown. Since the difference between the 
potentials is very small, panel (a), in panel (b) we show the tails of both potentials in the scale allowing to see the 
difference. Note that the existence of the potential (r), which is the phase equivalent to the bound state potential 

found in 0,0 with the same ANC, doesn't contradict to the inverse scattering theorem because the addition d ( ^| r - ) 
to the potential V N (r) asymptotically decays as exp(2 K a d,r) violating condition It doesn't allow one to use 
potential Vj (r) for analytical extrapolation of the scattering amplitude generated by this potential to the bound 
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FIG. 2: (Color online) The astrophysical factors 824(E) for the radiative capture a + d — t Li + 7. Black dots are data from 
3; black crosses are data from pj; black triangles are data from The red solid line is the 824(E) factor from 4|. The blue 
dotted line is S24(E1) factor, the blue dashed line is Sm(E2) and the blue solid line is our total astrophysical factor 824(E). 



state pole but we need this potential here only to generate the bound state wave function with correct amplitude of 
the tail, i.e. the ANC. 

Using the bound state wave functions generated by the potentials V N (r) and Vf(r) and the scattering a — d wave 
functions in the partial waves k — 1,2, we calculated the S 2 4(E) for the radiative capture a + d — > 6 Li + 7. The 
results of the calculations are shown in Fig. [2] The astrophysical factor with the bound state wave function generated 
by the potential V(r) is the 824(E) presented in [4J while our 824(E) is the one obtained using the bound state 
wave function generated by V\(r). At energies below the resonance practically the tails of the bound states wave 
functions do contribute to the matrix element. Since the square of the ANC in Q is higher than our one by 38%, 
correspondingly the astrophysical factor from [3| is systematically higher than our 824(E). Our calculations, definitely 
better reproduce the experimental data at energies larger than the resonance energies, where the calculations 
from [H clearly overestimate the data. At resonance energies both calculations reproduce the data @ very well. At 
astrophysically relevant energies E < 300 keV our 824(E) is lower than the one from Q by 38%. Finaly in Table Q] 
we present the a + d —> 6 Li + 7 reaction rates, which are also systematically lower than those presented in 

In this work we have demonstrated that a crucial quantity, which is necessary to pinpoint the 824(E) astrophysical 
factor, is the ANC for the virtual decay 6 Li — > a + d. Due to the peripheral character of the a + d — > 6 Li + 7 
direct radiative capture, this ANC determines the overall normalizatiion of the astrophysical factor at astrophysically 
relevant energies. From our calculations and Fig. [5] we can see that at low energies the contribution from the isospin 
forbidden El transition dominates over the allowed E2 transition. For example, at E = 70 keV, which is the most 
effective energy, the contribution from the El transition to the total 824(E) astrophysical factor is about 60%. Even 
at E ~ 100 keV the El transition contributes about 52% to the total astrophysical factor. Meantime, even if the 
Coulomb breakup of 6 Li would dominate, at E = 70 keV the El transition will be suppressed compared to the E2 by 
a factor of 60. It can hardly make possible to determine the total astrophysical factor from the 6 Li experiment. Since 
the ANC is the only crucial information needed to calculate the S 2 4(E) astrophysical factor at astrophysical energies, 
we call for more accurate measurements of the s-wave a — d elastic scattering phase shift at lower energies. It will 
help to extrapolate more accurately the data to the bound state pole to get more ANC for 6 Li — > a + d. Finally, we 
note that the problem of determination of the two-body bound state potential from the elastic scattering phase shift 
is quite important in different applications of nuclear reaction theory, in particular, in Faddeev approach for reactions 
with composite particles. 
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TABLE I: The rate for a + d — > 6 Li + 7 reaction calculated using our astrophysical S-factor for the temperature range 
10 6 K < T < 10 10 K. 



T 9 


N a {av) 


T 9 


N a (av) 




(cm 3 moV 1 s _1 ) 




(cm 3 mol -1 


0.001 


6.467 x 10 -30 


0.260 


6.823 x 10 -04 


0.002 


1.857 x 10~ 23 


0.270 


7.876 x 10~ 04 


0.003 


2.470 x 10~ 20 


0.280 


9.032 x 10~ 04 


0.004 


2.286 x 10~ 18 


0.290 


1.030 x 10~ 03 


0.005 


5.693 x 10~ 17 


0.300 


1.167 x 10 -03 


0.006 


6.592 x 10~ 16 


0.310 


1.317 x 10~° 3 


0.007 


4.651 x 10~ 15 


0.320 


1.478 x 10~° 3 


0.008 


2.327 x 10~ 14 


0.330 


1.652 x 10 -03 


0.009 


9.067 x 10~ 14 


0.340 


1.840 x 10~° 3 


0.010 


2.923 x 10~ 13 


0.350 


2.040 x 10~° 3 


0.011 


8.127 x 10~ 13 


0.360 


2.254 x 10 -03 


0.012 


2.008 x 10~ 12 


0.370 


2.482 x 10~° 3 


0.013 


4.508 x 10~ 12 


0.380 


2.725 x 10~ 03 


0.014 


9.343 x 10~ 12 


0.390 


2.983 x 10~ 03 


0.015 


1.811 x 10 -11 


0.400 


3.256 x 10~° 3 


0.016 


3.318 x 10" 11 


0.500 


6.930 x 10~ 03 


0.017 


5.787 x 10" 11 


0.600 


1.271 x 10~ 02 


0.018 


9.676 x lO" 11 


0.700 


2.148 x 10~° 2 


0.019 


1.559 x 10~ 1() 


0.800 


3.462 x 10~ 02 


0.020 


2.432 x 10~ 1() 


0.900 


5.385 x 10~ 02 


0.025 


1.538 x 10~ 09 


1.000 


8.079 x 10~° 2 


0.030 


6.277 x 10~ 09 


1.100 


1.166 x 10~ ()1 


0.035 


1.929 x 10~° 8 


1.200 


1.618 x 10~ 01 


0.040 


4.870 x 10~ 08 


1.300 


2.164 x 10~ 01 


0.045 


1.066 x 10~ 07 


1.400 


2.797 x 10~ ()1 


0.050 


2.093 x 10~° 7 


1.500 


3.508 x 10~ ()1 


0.060 


6.375 x 10~ ()7 


1.600 


4.288 x 10~ 01 


0.070 


1.554 x 10~° 6 


1.700 


5.126 x 10~ 01 


0.080 


3.245 x 10~ 06 


1.800 


6.002 x 10~ 01 


0.090 


6.057 x 10~° 6 


1.900 


6.915 x 10~ 01 


0.100 


1.038 x 10 -05 


2.000 


7.854 x 10~ 01 


0.110 


1.665 x 10" 08 


2.100 


8.808 x 10~ 01 


0.120 


2.533 x 10" os 


2.200 


9.773 x 10~ 01 


0.130 


3.690 x 10" os 


2.300 


1.074 x 10+ 00 


0.140 


5.185 x 10" 08 


2.400 


1.171 x 10+ 00 


0.150 


7.071 x 10" os 


2.500 


1.268 x 10+ 00 


0.160 


9.398 x 10" os 


3.000 


1.745 x 10+ 00 


0.170 


1.222 x 10~ 04 


3.500 


2.210 x 10+ 00 


0.180 


1.559 x 10~ 04 


4.000 


2.673 x 10+ 00 


0.190 


1.954 x 10~ 04 


4.500 


3.145 x 10+ 00 


0.200 


2.416 x 10~° 4 


5.000 


3.631 x 10+ 00 


0.210 


2.946 x 10~° 4 


6.000 


4.645 x 10+ 00 


0.220 


3.552 x 10~ 04 


7.000 


5.689 x 10+ 00 


0.230 


4.238 x 10~ 04 


8.000 


6.725 x 10 +0<) 


0.240 


5.008 x 10~ 04 


9.000 


7.723 x 10+ 00 


0.250 


5.868 x 10~ 04 


10.00 


8.664 x 10+ 00 
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